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Abstract. We present preliminary results of our numerical study of the critical dynamics of perco- 
lation observables for the two-dimensional Ising model. We consider the (Monte-Carlo) short-time 
evolution of the system obtained with a local heat-bath method and with the global Swendsen-Wang 
algorithm. In both cases, we find qualitatively different dynamic behaviors for the magnetization 
and £2, the order parameter of the percolation transition. This may have implications for the recent 
attempts to describe the dynamics of the QCD phase transition using cluster observables. 



The study of the dynamic critical behavior of simple statistical models might be of rel- 
evance for understanding non-equilibrium effects in hot QCD, such as the effects due 
to heating and cooling of matter produced in heavy-ion collisions. The possible connec- 
tion yj] between the deconfinement transition in QCD and the percolation phenomenon 
J2] has received renewed attention recently J3] and the dynamics of cluster observables 
has been investigated using hysteresis methods [4]. As a preparation for the study of the 
continuous-spin 0(A) model, whose magnetic transition is expected to be in the same 
universality class of the chiral phase transition in two-flavor QCD, we consider here the 
short-time dynamics of the 2d Ising model and focus on the dynamic critical behavior 
of percolation observables. 

For many physical systems, a suitable definition of cluster provides a mapping of the 
physical phase transition into the geometric problem of percolation, allowing a better 
understanding of how the transition is induced in the system. For Ising and O(N) spin 
models this mapping is well understood 0,1a], whereas in QCD it may be harder to 
define, even in the pure-gauge case [3]. Given a definition for a cluster on the lattice, the 
order parameter in percolation theory is the stress of the percolating cluster Q., defined 



where A is the volume of the percolating cluster and T c is the critical temperature. 



Note that a percolating cluster is a set of spins connected from the first to the last row of the lattice. 



INTRODUCTION 



by 




(i) 



SHORT-TIME (MONTE CARLO) DYNAMICS 



Using renormalization-group theory, it can be shown [7] that the early time evolution 
of an order parameter (e.g. the magnetization M) already displays universal critical 
behavior, given by 

M(t,£,m ,L) =L-P lv Jt{tL-\£l}l v ,m Q L Xa ), (2) 

where £ = (T — T c )/T c and ^ is a universal function. We thus expect for T = T c 
a power-law behavior at early times M(?) e _>o ~ mQt e , with = (xq — j8/v)/z. In 
principle, we would assume that the percolation order parameter £1 should have a similar 
behavior. The time evolution for the heat-bath and cluster algorithms is described below. 

The heat-bath dynamics |0] consists in choosing the two possible directions of each 
Ising spin according to the exact conditional probability given by its nearest neighbors. 
Each spin S,- is chosen UP with probability pi , or DOWN with probability 1 — pu where 



" l+exp(-2/3E ; 5 7 ) 

and the sum is over nearest neighbors of Sj. After a certain number of iterations the 
spin configuration obeys the Boltzmann distribution. In the heat-bath method, since the 
updates are local, this transient time becomes considerably large at criticality. 

The Swendsen-Wang (cluster) algorithm is obtained from the Ising-model Hamilto- 
nian by writing the partition function as 

C» y =i 1 r«, 7 =o 1 

^=ii n^-pij) > w 

{S}{n} { J I (ij) J 

where py = 1 — exp(— 27/3) is the probability of having a link between two nearest- 
neighbor sites of equal spin value. This link is represented by nij and determines the 
clusters that will be associated with percolation at the critical temperature [9]. The 
dynamics then consists of global moves in which the spins of a cluster are flipped 
together, with probability 1/2. 



NUMERICAL RESULTS 

In order to study the short-time dynamics we simulate at T = T c and force the system to 
have an initial magnetization mq. We let the system evolve in time and look for power- 
law behavior of the order parameters M and Q. as a function of the (Monte Carlo) time. 
Each temporal sequence is generated from a different random seed, i.e. each sequence 
has a different initial spin configuration. The time history is then obtained from an 
average over all the generated sequences. We have performed Monte Carlo simulations 
with 50000 seeds and 5000 sweeps for three initial magnetizations (mo = 0.02, 0.03, 
0.04), for six lattice volumes, for both heat-bath and Swendsen-Wang algorithms. 
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Heat Bath, L=100, 50000 seeds, 500 sweeps (Monte Carlo time) 




FIGURE 1. Plot of the early time evolution of the magnetic (M) and percolation (£1) order parameters 
for the heat-bath case. 



We obtain that a power-law fit works very well for M, yielding the literature value ifioll 
for the exponent 6. However, as can be seen in Fig.[H the percolation order parameter £1 
does not show a power-law behavior, being consistent with an exponential behavior of 
the type exp(— v/t), where x is essentially the relaxation time to equilibrium. This seems 
to hold also for the case of the Swendsen-Wang dynamics. We thus see that although 
the equilibrium behaviors of M and Q. are equivalent, the two order parameters show 
qualitatively different dynamic critical behavior. 
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